Dynamics of elastic rods in perfect friction contact 
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One of the most challenging and basic problems in elastic rod dynamics is a description of rods 
in contact that prevents any unphysical self-intersections. Most previous works addressed this issue 
through the introduction of short-range potentials. We study the dynamics of elastic rods with 
perfect rolling contact which is physically relevant for rods with rough surface, and cannot be 
described by any kind of potential. We derive the equations of motion and show that the system 
is essentially non-linear due to the moving contact position, resulting in a surprisingly complex 
behavior of the system. 
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Introduction Take two rubber strings, stretch them a 
bit and cross them so they remain in contact, as shown 
on Figure 1. As long as the deformations are not too 
large, the strings will roll at the contact without sliding. 
Of course, that simple experiment is dominated by the 
energy loss from the internal deformation of strings at the 
contact point; however, improving the quality of strings 
will allow them to oscillate for a reasonable time before 
the energy loss takes over. Interestingly, this simple and 
familiar experiment has deep mathematical and physical 
implications that go beyond a toy problem. 

There are many objects that can be represented as 
long elastic rods, from a rubber hose to DNA molecules. 
Typically, if these objects are put in a confined space, or 
undergo other non-trivial dynamics, self-contact of these 
rods typically appears. The true dynamics will not al- 
low the rod to pass through itself, and it must preserve 
the side of contact under dynamics. While it is gener- 
ally accepted that something like DNA at contact will 
slide freely, the dynamics of other molecules like den- 
dronized polymers (DP) may be different. These com- 
pound molecular structures are formed by assembling 
multiple dendrons that are each connected by its base 
to a long polymeric backbone [1]. A simplified, coarse- 




FIG. 1: A sketch of two strings in contact. The centers of 
disks in contact are marked r\ and T2, the line connecting 
them passes through the contact point. Reaction force of the 
contact is also shown as A. 



grained rod model of such polymers must take into ac- 
count the rough surface formed by tree-like structures 
that is likely to generate perfect rolling contact. 

Efficient numerical methods have been developed re- 
cently to deal with the self-contact forces of rods using 
short-range repulsive potential for statics [ ] and dynam- 
ics [3-5]. Another avenue of studies of stationary states 
of elastic rods with self-contact [6-10] explicitly computes 
the contact forces from the existence of constraints. In 
our case, the rolling contact comes from friction, which 
does not admit any potential description. The extension 
of the latter approach to include the dynamics, and es- 
pecially the rolling constraint is difficult. We also note 
that the true motion may be a combination of rolling and 
sliding friction, but in the absence of a consistent theory 
for rolling slippage we shall concentrate on the perfect 
rolling only. Our solution of the problem of rolling con- 
tact is similar in spirit of the second approach, as the 
contact force comes naturally from the constraint and 
not from short-range forces. 

Setup of the problem As is well known [ ], the prob- 
lem of rolling motion is essentially non-holonomic and, in 
general, cannot be represented from the potential point 
of view. We shall note that it is possible to define poten- 
tial approach to some non-holonomic systems [ ] . How- 
ever, these cases seem to be more of an exception rather 
than a rule. For the problem in hand, we proceed with 
the Lagrange-d'Alembert (LdA) principle which is the 
fundamental method for the treatment of non-holonomic 
systems [11, 13]. Similar approach has been recently used 
to describe the motion of an elastic rod rolling on a plane 
[14]. In order to utilize LdA, we recast the motion of the 
rods in variational setting through the Simo-Marsden- 
Krishnaprasad (SMK) rod theory [15], which allows a 
variational formulation of string dynamics [ ]. The rod 
is parameterized by a coordinate s that does not have to 
be arc length. We fix a reference frame and measure the 
local position r(s, t) and orientation A(s, t) at time t with 
respect to the fixed frame. SMK theory formulates the 
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dynamics in terms of the variables that do not depend 
on the choice of basic frame: 

T = A-V , n = A" 1 A' , 7 = A- X r , w = A" 1 A . (1) 

Here, the prime denotes the partial derivative with re- 
spect to s and dot the derivative with respect to t. The 
physical meaning of 7 and uj is the linear and angular 
velocity in the body frame of reference, and T and ft 
the corresponding deformation rate. If £(r,f2,7,u;) is 
the Lagrangian, then the equations of motion for a free 
elastic rod are given by the variational principle [16] 



S J £(T,n,<y,u;)dsdt = 0, 



(2) 



appropriately computing the variations of variables 
(T,f2,7,u;). For purely elastic rods, £ is a quadratic 
function; more complicated expressions for £ are also pos- 
sible. In the case of constraints, this variational principle 
(2) has to be modified according to the LdA approach as 
follows. Assume, for simplicity, that the undeformed rod 
has a circular cross section and the acting forces are small 
enough so that the cross-section remains circular even at 
contact. Thus, if the strands at contact can be approx- 
imated by touching circular cylinders, then the contact 
point is always located at c(t) = (ri + r 2 )/2, and the vec- 
tor from the center of the cylinder Vi(t) to the contact 
point is ±(r 2 — ri)/2. Here, the index i = 1,2 means 
evaluation at s = Si marking the disks at contact. Since 
the angular velocity in the fixed frame is A^A" 1 , i = 1, 2, 
then the velocity of the material point associated with 
contact, also in the fixed frame, is 

ri + ^ Ai AJ- 1 (r 2 - ri) = f 2 - ^ A 2 A^ 1 (r 2 - ri) (3) 

Note that it is not the same velocity as the velocity of the 
contact point, which will include the full derivatives, and, 
as a consequence, S{. This condition can be reformulated 
in terms of invariant variables by multiplying (3) by A^ 1 : 

7i + x «12 = £1272 - \ (£12^2) x k 12 , (4) 

where £i 2 := Aj~ 1 A2 is the relative orientation, K12 = 
A^ x (r 2 — Ti), and other invariant variables are defined 
as in (4), with the index i meaning evaluation at s = s$. 
Note that due to the uniformity of the strand, and the 
assumption of circular cross-section, the Lagrangian does 
not depend explicitly on Si and S{. Then, LdA principle 
states that one has to replace time derivatives only in (4) 
by (5-variations and use that expression as an additional 
constraint on (ST, Sfl, £7, Suj) in (2). For example, 

7i = Ar'fi A^Sr(s)S Sl , (5) 

and similarly for other variables. Here we have denoted, 
for shortness, S Si := S(s — Si). There is an unfortunate 
collision of notation in S between the variational deriva- 
tives and Dirac's (5-functions; in this paper, the (5-function 
always has a subscript. Note that one does not enforce 



the condition (4) as it stands: this approach leads to the 
so-called vakonomic approach ['. ] which can be shown 
to give incorrect equations of motion for this case. 

Denote by X(t) a vector that enforces the LdA con- 
straints, and by D/Dt := d t -f u>x, D/Ds := d s +fJx the 
full t- and s- derivatives. We get the following equations 
for the motion of strings with rolling contact: 



D d£ 
Dt£hy 
D d£ 



D d£ 
DsdT 

D S£ 



xs s 



■ ^ 12 1 \5 S2 
d£ 



^ X7+ ar xr 



i/«i2 x X5 Sl + ^2 (k 12 x \j5 S2 . 



(6) 



(7) 



The physical meaning of A as the force due to the con- 
straint is now evident from the linear momentum (6). 
Correspondingly, the A terms in equation (7) are iden- 
tified as torques acting on the contact point due to the 
presence of the constraint. These equations have to be 
augmented by the compatibility conditions 

f = 7'-u>xr + f*X7 

(8) 



ft 



x n. 



We still have to close the system by computing the equa- 
tions of motion for the contact points s^, which is done 
from the tangency conditions stating that the strand at 
a contact point, which is locally a cylinder, touches itself 
tangentially and there is no intersection of the material. 
These conditions state that 



Ki2(si,s 2 ,t) ■ ri(si,t) = 0, 
«i2(si, s 2 ,i) • & 2 r 2 (s 2 , t) = 0. 



(9) 



Note that these expressions do not contain any time 
derivatives of the variables and are thus holonomic. They 
also imply that the distance between the centers of disks 
in contact \r± — r 2 | is preserved. Since the Lagrangian £ 
does not depend on and s^, they can be imposed after 
the equations have been derived. In principle, such con- 
ditions already determine S{ through an implicit relation; 
however, such relation is difficult to use. Instead, we take 
the time derivatives of the tangency conditions to obtain 
a relationship 

A-(si, 5 2 ) T + v = 0, (10) 

where the 2x2 matrix A and 2- vector v depends on the 
dynamical properties at contact, with det(A) =^ when 
the rods are not locally parallel at the contact. 
Discrete strands in contact It is also interesting to 
consider the application of this theory to discrete, chain- 
like elastic structures in contact. In that case, we need 
to clarify the physical meaning of the (^-function at the 
contact position. Apart from its physical relevance, this 
consideration is also useful for consistent numerical dis- 
cretization of (6,7). Here, care must be taken in deriv- 
ing the equations of motion without breaking their vari- 
ational structure [18]. Suppose that we have a string 
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consisting of discrete set of points along the line, s = 
with k being integer. If the orientation and position of 
a material frame at s = s k are given by an orientation 
matrix and a vector r/~, the invariant variables are 
p k = Aj^Afc+i and q k = Aj^rfc+i - r k ). The purely 
elastic Lagrangian is £ = £(uJk,7kiPkiQk)' One also 
needs to define a "smeared-out" version of (4): 

a k(l k + ^fimUk X K km ) 

-OLkPm{ikmlm ^(6cm^m) X K km ) 0, (11) 

(summation over fc,m). Here, we have defined the aver- 
aging coefficients: a k '•= G(si — s k ) /3 k '•= G(s2 — s k ) 
arising from a "bump" function G(s) that rapidly de- 
cays away from s = 0, and £ km := A" 1 A m , n krn := 
A^ 1 (r m — r/c). Note that the positions of the contact 
is defined by the continuous variables s\ and 52. The 
physical meaning of (11) consists in spreading the point 
wise contact condition (4) to a few neighboring points 
surrounding the contact. Then, the LdA principle gives 
a discrete analogue of (6,7): 

D di , d£ d£ n , w . 

mwr^-^+w^ (12) 

D di d£ _! _! ^ _ a/ 

\ (afcfc^ m xA+fta m ^(K mfc xA)) . (13) 

It is interesting to note that for a Lagrangian that is 
quadratic in all variables (e.g., linear elasticity), (12,13) 
allow explicit calculation of the constraint A. These 
equations have to be augmented by conditions for the 
variables si and s 2 . The discrete version of (9) is 

Yjhn a lPm^lm ■ T[ = and Yjhn a lPm^lm ' (irnXrn = 0. 

Differentiating these expressions with respect to t, we 
again obtain the equations governing the evolution of 
derivatives si and S2 as A^- (si , S2) T + Vd — 0, where the 
matrix and the vector depend on the dynamical 
variables close to the contact. 

Linear strings in contact One may wonder if the 
equations of motion we have derived allow to deduce 
analytical expressions for the propagation of the dis- 
turbances along the rods at contact. The answer to 
this question is, unfortunately, no: the contact condi- 
tion makes the disturbances essentially nonlinear. Let us 
consider two strings in contact, and denote for shortness 
a = (7,u>) T and A = (r,f2) T . For linear elastic mate- 
rials d£/daL = Va and 81/ d A = -QA, where V and Q 
are 6x6 matrices. The linearized compatibility conditions 
(8) allow to introduce a vector potential as a = 9^0, 
A = d s <j). Neglecting all nonlinear terms in the dynamic 
variables and assuming that the rod is naturally straight 
in its undeformed state, we can transform (6,7) into a 
vector wave equation [19] 



B 2 ct) B 2 (t> ( Id \ ~ , t J -Id \ . . nA . 
The condition (4) can be expressed in this vector form as 

(ld,-iKi2x) T t (5i,t)=(^i2,iKi2X^i 2 ) T ^(s2,t) (15) 

Thus, the evolution of small disturbances on the rod 
is governed by linear equations (14), linear rolling con- 
straint (15) and nonlinear evolution equations for si^ 
(10). We can illustrate the complexity of this problem 
on a pedagogical simplistic example of two straight rods 
in contact with only one mode being relevant in (14) for 
each rod. Let us call these modes u(x, t) for the first rod 
and v(y,t) for the second rod, respectively. Such a sim- 
ple model illustrates the full complexity of the problem 
without the unnecessary complications of higher dimen- 
sionality of <p. Physically, such a mode can be realized 
for a rod with special elastic and inertia matrices V and 
Q, e.g. for rods made out of composite material. In this 
case, the rolling constraint and the motion of the contact 
points are written simply as 

ut(si) = Fv t (s 2 ) , ^ = GiUt(si) , i = 1, 2 , (16) 
where F and Gi are constants depending on the material 
parameters and the base state of the rods. The equations 
of motion for (u, v) in this reduced setting are: 

u t t ~ c 2 u xx + X5 8l = , v t t - c 2 Vyy - XFS S2 = 0, (17) 
where A enforces the first constraint of (16). 

A discrete version of these equations can be derived, 
similarly to the full equation described above in (13). We 
shall emphasize that the complexity caused by the non- 
linear contact conditions is the same in the full equation 
(14) and its one-dimensional counterpart (17). While the 
complex dynamics caused by the nonlinear rod equations 
have been well studied, as far as we are aware, there has 
been no work on the complex dynamics caused by the 
rolling contact condition. In the absence of constraints, 
the wave equations provide a simple harmonic oscillation 
of the string. However, when the constraint is present, 
the motion of the string is challenging and complex. 
Wave equation: exact solution and contact chaos 
In the case when u(x, t) and v(y, t) have periodic bound- 
ary conditions on an interval of the same length, an ex- 
act solution of the equation (17) can be found, in spite 
of the contact nonlinearity. Let us consider the case 
when the strings are identical, so G\ = G2 = G and 
F = 1. In that particular case, from (16), we conclude 
that 52 — si = d =const. By shifting the coordinate y in 
the ^-equation of (17) as y = x+d, using periodic bound- 
ary conditions and denoting the new coordinate x so both 
u and v depend on the same variables, we can choose 
si(0) s 2 (0). From si = s 2 , we get si(t) s 2 (t) = s(t). 
Define P — u + v and M = u — v, satisfying 

Ptt ~ c 2 P xx = , Mtt ~ c 2 M xx = 2X5 (x - q(t)) (18) 
M t (S) = 0, s= l -GP t (s). (19) 
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As we see, equation (18) for P(x, t) is now decoupled and 
can be solved independently from the other equations as 
a simple wave equation on the circle. Then, s(t) is ob- 
tained by solving the first-order differential equation (19). 
To find the solutions u(x,t) and v(x,t), one has to solve 
equation (18) with prescribed s(t) which is a much sim- 
pler problem. The solution P{x,t) can be any function 
satisfying a wave equation with periodic boundary con- 
ditions. The equation for s(t) is then s = GP(s(t),t)/2 
which is a nonlinear, non- autonomous equation; it is well 
known that such equations lead to rather complex behav- 
ior for s(t) even if P(x,t) is a simple function. However, 
this behavior is not truly chaotic in our opinion. 

The situation is much more interesting for the case 
of more realistic fixed boundary conditions for u and v. 
This would be the case of two rods with fixed ends in 
contact. Such motion, as one can show from equations 
(17), conserves energy; however, in reality, friction with 
air and, more importantly, rolling friction will lead to 
energy dissipation. It is nevertheless interesting to see 
the structure of that dynamics, with the nonlinearity ob- 
tained only from the contact condition. As we see from 
Figure 2, the system produces a complex spatio-temporal 
dynamics of both strings. It is also relevant to present 
another measure of complexity, computed from the dy- 
namics of the base harmonic of u(x,t), call it ui(t), as a 
function of t. If the string were vibrating in the air, the 
sound sufficiently far away from the string will primarily 
contain the contribution from the first harmonic. On Fig- 
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FIG. 2: Left: the spatio-temporal evolution of one of the 
strings in contact, with the red line marking the point of con- 
tact. Right: spectrum S(cj) of time signal produced by the 
lead harmonic (in x) of u(x,t). 

ure 2, right, we plot the time spectrum S(oo) as a function 
of temporal frequent lj, obtained from the time signal of 
the first harmonic u\(t). Starting with u(x, 0) = sinx, a 
linear rod with < x < 2ir without rolling contact will 
generate a purely monochromatic sound; however, when 
the contact is present, there is a persistence of high over- 
tones to the signal. We have found out that the chaotic 
behavior caused by the contact condition persists for all 
initial conditions we have tried. We believe that the ap- 



pearance of the chaos due to the contact is highly inter- 
esting and important for many physical applications, and 
yet it has not been discussed in previous literature. 
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